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The oblique reflection of VLF radio waves from a continuously stratified ionized medium 
is considered. The profile of the effective conductivity is taken to be of an exponential 
form. This is a fair representation for the actual D layer of the ionosphere. It is shown 
that the gradient of the conductivity change has a marked effect on the reflection 
characteristics. 



1. Introduction 

In the study of VLF radio wave propagation it is 
often assumed that the ionosphere can be regarded 
as a sharply bounded medium on its underside. 
This step may be considered reasonable in view of 
the long wavelength and the relatively rapid change 
of the electron density of the D region. Further- 
more, experimental data are often in accord with 
calculations based on this model. The general 
agreement is particularly good at highly oblique 
incidence provided that appropriate corrections are 
made for earth curvature. Nevertheless, there are 
many occasions when the sharply bounded model 
appears to be inadequate. Therefore, it is worth- 
while to consider a more realistic model of the lower 
ionosphere. 

On examining much of the recent literature on the 
characteristics of the lower D layer, it appears that 
the effective dielecti^ constant of the medium can 
be well approximated by an exponential function. 
Then, to within this approximation, the relative 
permittivity may be written in the form 

K(z)=K (l-ij- exp jsA (1) 

where K is a reference permittivity, and L and /3 
are constants. The level 2=0 may be denned as the 
reference level and thus 



which is a familiar form. 



(2) 
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Under the assumptions that the angular frequency 
o) is much less than both the collision frequency v 
and the plasma frequency co , it is known that [Wait, 
1962] 



Ko=l and L=— > 
Wo 



(3) 



provided that the earth's magnetic field can be 
neglected. The latter assumption is strictly valid 
only if v is somewhat greater than the sz;yrofrequency 
co//. In actual cases, however, the isotropic assump- 
tion is useful even when v is of the order of co#, 
provided the magnetic field is not transverse to the 
direction of propagation [Crombie, 1961; Jollier and 
Harper, 19621. For the purposes of the present 
paper, the influence of the terrestrial magnetic field 
is not considered. 

The constant (3 in the exponent of the equation for 
K(z) is a measure of the sharpness of the gradient. 
For example, when /3— 1 km -1 , it means that the 
ratio o)l I V or Njv increases by a factor 2.71 for each 
km of vertical height. The best available informa- 
tion on expected values of p can be found in the 
recent work on pulse cross modulation [Barrington 
et al., 1962] and from recent rocket measurements 
[Kane, 1962]. From these results it appears that 
P is of the order of % km" 1 for quiet daytime 
conditions although it may differ by a factor of two 
or more at certain times. In the present study it 
appears to be desirable to allow to vary from 0.3 
to 3.0 to encompass most cases of interest. 

2. Background Theory 

The propagation of electromagnetic waves in a 
medium whose permittivity e(z) varies in an expo- 
nential manner has been studied extensively. The 
usual case considered is for horizontal polarization 
where the electric vector is parallel to the stratifica- 
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tion. For example, if the electric vector has only a 
y component E v , it is a simple matter to show that 
[Wait, 1962] 



\Jz 



;+kKK(g)-S*) \E y =0, 



] 



(4) 



when the field varies in the x direction according to 
exp (—ikSx), where S is a dimensionless constant. 
Here fc=V € oMo^ is the free-space wave number and 
e and /z c are the permittivity and permeability of 
free space, respectively. 

At a sufficiently large negative value of z, the 
relative permittivity K(z) becomes unity, and thus, 



E y satisfies 



(£+ w ) 



E„=0, 



(5) 



where C 2 = l—S' 2 . The general solution of this 
equation has the form 



E y =E {e- tkCz + R h e+ tkCz )e- 



(6) 



when E and R h 



^0 emu. lih are constants. Recognizing that 
the solution of (4) must reduce to (5), for z tending 
to — co ? it may be shown that [Wait, 1962] 



«»= 



~(ITr- 



1 



-ay* ( Vo y. 



(7) 



where v Q =2ikC/p. 
The term Etf- ikCz '< 



may be regarded as an 



incident wave whose direction of propagation makes 
angle 6 with the z axis where C=cos 6 or S=sin 6. 
Then E Re +ikCz e~" ikSx can be interpreted as a reflected 
wave and R can be defined as the reflection coeffi- 
cient of the exponential layer. It should be stressed 
that this reflection coefficient, while referred to the 
level 2=0, is only valid in the free space region at 
large negative values of z. 

The reflection coefficient R h , given by (7), may be 
conveniently written in the form 



where 



<*>=- 



U 



R h =exp (-g tf) exp **, (8) 

^-^>+ 2 -'[(-^>} 



where \=2irjk is the free-space wavelength. Thus, 
the amplitude of the reflection coefficient is given by 
the very simple form 



i*.i=«*(-g* 



(9) 



which is independent of L. On the other hand, the 
phase factor <£> is relatively complicated since it in- 
volves the factorial function of imaginary argument. 
Nevertheless, (8), for the complex value of R h can 



be used to obtain quantitative results for exponential- 
type layers. 

3. Extension to Vertical Polarization 

In the VLF radio problem one is mainly interested 
in vertical polarization such that the magnetic vector 
is parallel to the stratification. For example, if the 
magnetic field lias only a y component, H VJ and 
assuming again that the fields vary in the x direction 
according to exp (—ikSx), it follows that [Wait, 1962] 

Intrinsically, this is a more complicated form than 
(4) for horizontal polarization. 

In the region of large negative z the (10) reduces 
to the elementary form 



(jp+Fc)ff„=0, (11) 



and thus 



fl;=flr (e-«»«+B f «+«»«)«~ ttte , 



(12) 



where H G is a constant and R v is, by definition, the 
reflection coefficient for vertical polarization. 

Unfortunately, for the exponentially varying 
permittivity of the form defined by (1), it does not 
appear to be possible to obtain a closed form expres- 
sion for R v . For this reason, the required quantita- 
tive results were obtained by a numerical method. 
Essentially, the procedure is to replace the contin- 
uous K(z) profile by a finite number of steps. The 
situation is illustrated in figure 1 where the function 



lfe) = Z exp ^' 



(13) 



is shown plotted versus z along with its step approxi- 
mation . Thus , between the limits z = — z and z = -f- T, 
the medium is divided into M homogeneous layers of 



width hi y h 2 



hx 



h M . The value of 



L(z) in each of these slabs is then replaced by a 




Figure 1. The step approximation to an exponential profile. 
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constant value L m . The method is really equivalent 
to the usual method [Budden, 1961] of numerically 
integrating the basic differential equations. 

The problem now boils down to finding the wave 
impedance Z } at z= — z in terms of the properties 

W 
of all the individual slabs. Here Z 1 = T j which is to 

be evaluated at z— — z . The reflection coefficient R v , 
which is referred to the level 2=0, is then found 
from 



R r 



C-ZJ Vo 

'C+ZJvo 



exp 



[-SH 



(14) 



for z= — z , and where rj =-\ , JjL /e = 120 w. Here z is 
chosen sufficiently large that 1/L(z) may be regarded 
as zero. (Strictly speaking, z ^oo.) 

4. Iterative Process 

From the theory of wave propagation in stratified 
media [e.g., Wait, 1962], it is known that Z, can be 
obtained by a series of iterative processes. Thus 



where 



Qr- 



Z 1 =K 1 Q l , 

l-\-K 2ii Q 2 tanh D x 



n K 3 , 2 Q 3 +Un\i D 2 
V2 l+K s , 2 Q 3 t&nhD 2 



Vm l+K m+Um Q m+1 t&nh D m ' 



(15) 
(16) 

(17) 

(18) 



and so on. The various factors arc defined below. 
The parameter D m which is a measure of the thick- 
ness of an individual slab is defined by 



D„. 



i'Zwh,, 



N r , 



S 2 \x 



0-0 



(19) 



where N m is the refractive index of the mth slab 
given by 2 



N m = 



(-0* 



The quantity K m +i,m is the ratio of the wave im- 
pedances of layer m-j-1 and layer m. Thus 



J7- -A-m+l 



(20) 



where 



K„, 






S 2 \» 



(-i) 



2 The square roots are defined such that D m = ( T jN m if AV-* 00 , while 
D m =(^~ n ^\C if 2V M ->1. Furthermore, N m ^0/L m )^ 2 e~ iT/4 it \N m \»l and 



< >.o 
A T m = lif L m -> 



is the wave impedance in the mth layer. Alternately, 



■K-m + l.r. 



('-£) 
"0-i4 



C 2 - 



Arc + 1 

rn !L 

J 



(21) 



Finally, because the quantities Q m are given in 
terms of Q m +\ it is necessary to know the initial value 
Qm+i at z=T. For example, 



Qm- 



Km+i, mQm+i+ -tanh D M 
~l+K M +i,M[tnnh D m ]Qm+i 



(22) 



and similarly, Q M -\ is expressed in terms of Q M . The 
process is continued until Q x is obtained. 

In some problems of this kind it is not necessary 
to know precisely the value of Q M +\- However, the 
economy of the calculations is greatly improved if a 
good starting value of Q M +i is known. Recognizing 
that the wave impedance at z=T is given by 



y 



K 



M+l 



1)0 



& 



(23) 



and utilizing the condition |Af M+ i|>>l, it follows 
from previous work [Wait, 1962J that 



K 



Qm+i* 



K x 



/- 2t Ar 2 



V* T" N M +i 

An 



(24) 



where Kq and K Y are modified Bessel functions. 

In a practical sequence of calculations it is neces- 
sary to choose T large enough that the final result for 
R v is insensitive to further changes in T. For ex- 
ample, if ]S=1 km -1 and X =15 km, it was found that 
7=4 km was sufficiently large to obtain four-figure 
accuracy in R v . Furthermore, for some values of f$ 
and Xo, it was found that R v approaches a constant as 
zo was increased to 20 km. In general, for smaller 
values of /3 it was necessary to increase T and Zq to 
larger values in order to achieve stability of the final 
results. The width of the steps, h m , must also be 
chosen sufficiently small to achieve an adequate 
simulation of the smooth profile. Generally, it was 
found that \D m \<^lO~ 2 was a satisfactory criterion. 

5. Discussion of Reflection Coefficient 
Calculations 

Using the basic definition of the reflection coeffi- 
cient R v , given by (14), the amplitude and phase of 
R v were calculated for a wide range of conditions. 
While, strictly speaking, R v is defined as the limit 
for z= — Zo, where zq—>™, it is of interest to first 
demonstrate how R D approaches this limit. An 
example of such calculations is shown in figures 2a 
and 2b. Here (7=0.1, corresponding to highly 
oblique incidence, and X =15 km, or /=20 kc/s. 
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Figure 2. The reflection coefficient as a function of height for 
an exponential profile and a fixed angle of incidence. 



Various values of are shown on the curves. It is 
immediately apparent that for sufficiently negative 
values of z (i.e., sufficiently far below the layer), 
both the amplitude and phase of R v approach a 
limit. 

Actually, the curves in figures 2a and 2b have more 
than just a mathematical interest. Physically, the 
R v , corresponding to z= — 2o where z is finite, is the 
reflection coefficient for a permittivity profile of the 
form 



K(z)=K [l-~e^lior 
=0 for z<C — z . 



O — 2o 



(25) 



Therefore, these curves describe reflection from a 
sharply bounded ionosphere which behaves expo- 
nentially above its lower edge. Of course, as zo 
tends to oo ? the discontinuous profile becomes a 
continuous exponential. It is rather interesting to 
observe that as the discontinuity is moved from 
above the reference level (2=0) to below, the 
reflection coefficient passes through a minimum. 
This is related to a Brewster angle phenomenon 
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Figure 3. The amplitude of the reflection coefficients for 
vertical and horizontal polarization as a function of C which 
is the cosine of the angle of incidence. 



since to the right of this minimum in figure 2a the 
reflection is metalliclike while, to the left, the 
reflection is dielectriclike. 

For the remainder of the present paper, attention 
will be confined to the limiting value of R v far below 
the reference level 2=0 (i.e., 2^ — 2 =~°°)- On 
the other hand, it should be emphasized that R v is 
always referred to the level z=0. 

The magnitude of the reflection coefficients R v and 
R n are shown plotted in figure 3 for vertical and 
horizontal polarization, respectively. The values 
of \B V \ were obtained from (14) using the multislab 
model described. The values of \R h \ were calculated 
directly from (9) and, on this log-linear scale they 
are merely straight lines. As an important check, 
the value of \R h \ was also obtained from a multislab 
model using the same procedure as described for 
vertical polarization. Within four-figure accuracy, 
the values of \R h \ obtained by the two methods were 
the same. 
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A number of significant features are evident in 
figure 3. In the first place, |B r | exhibits a Brewster 
angle phenomenon, provided p is sufficiently large. 
Here the reflection process is dielectriclike at grazing 
angles and is more or less metalliclike at normal 
incidence. Of course, \R n \ does not exhibit this 
phenomenon. However, for very small values of ft, 
corresponding to a relatively slowly varying medium, 
the curves for \R h \ and \R P \ become rather close to 
one another. This is consistent with the optical 
behavior of waves in an inhomogeneous medium. 
In fact, by a direct application of the phase integral 
method of Eckerskey, the equation for \R h \ is found 
also to be applicable to \R V \. The use of the phase 
integral method for such applications is necessarily 
restricted to slowly varying media [Wait, 1962]. It 
is apparent, that for VLF radio waves, where p is of 
the order of 0.5, the phase integral method is 
inapplicable to vertical polarization. 

Another interesting feature of figure 3 is the near 
linear dependence of all the curves for small values 
of C. Fortunately, it is just these values which are 
important in the long-distance propagation of VLF 
radio waves. The linearity of the phase curves for 
Rj, are indicated in figure 4 for the same conditions. 
Here, it is evident that they all approach —180° at 
grazing incidence. The phase curves for R^ exhibit 
a similar property but they are not shown here since 
they have only an academic interest at VLF. 

The variation of the reflection coefficient R„ as a 
function of the gradient parameter p is illustrated 
in figures 5a and 5b at oblique incidence. It is 
rather remarkable that \R V \ is relatively insensitive 
to P if it is in the range from 0.7 to 3.0. Further- 
more, it appears that \R V \ has a broad maximum for 
P approximately equal to 1.2. 

The general behavior of the amplitude and phase 
of the reflection coefficients at highly oblique inci- 
dence suggests that, if R v is written in the form 

R v =— exp(aC), 
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the function a should be almost independent of C. 



Writing 



a=ai-\-ia 2f 



where a x and a 2 are real, it is a simple matter to 
compute the complex coefficient a from the numerical 
data of R v . The results are shown in figures 6a and 
6b where ( — a x ) and a 2 , respectively, are plotted as 
a function of C in the range 0.05 to 0.30. It is 
apparent that over this range of C, the coefficient a 
can be regarded essentially as a complex constant. 
For all the results given in the foregoing figures, it 
has been assumed that L=l/2 and A = 15. Actually, 
the results for the magnitude of R v and R h do not 
depend on L. In fact, \R V \ and \R h \ are identical for 
/3\o= constant. Furthermore, the phase of the 
reflection coefficients is also simply related although 
the situation is slightly complicated by the choice of 
the reference level where L is required to have a 
special value. In fact, it is convenient to choose the 
reference level (i.e., 2=0) so that 



L=- 



co 



"2 X X„ 



In this way, the results can be readily compared as a 
function of frequency or wavelength. At the refer- 
ence level z=0, the effective conductivity parameter 
w, or u>l/v has the value 2.51 X10 5 . 
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Figure 5. The reflection coefficient as a function of the 
gradient parameter /3 for an exponential profile. 
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defined by R v =—exp(aC) for X =15 and various values of /3. 



For wavelengths other than X =15, it is also found 
that R v may be approximated by the function— 
exp(aC) where a is approximately a complex con- 
stant. To illustrate the wavelength dependence, 
— «i and a 2 are plotted as a function of X in figures 
7a and 7b, respectively. For these curves, C is 
chosen to be 0.16; the corresponding curves for other 
values of C in the range 0.10 to 0.20 are almost 
indistinguishable. The ordinate in figure 7a is 
simply related to the magnitude of the reflection 
coefficient and, thus, small values of — a x are associ- 
ated with high reflection coefficients. It is apparent 
that for small values of fi corresponding to a diffuse 
layer, the reflection coefficient becomes very small for 
the shorter wavelength. On the other hand, for a 



rapidly varying layer, corresponding to large values of 
0, the reflection coefficient decreases with increasing 
wavelength. For intermediate values of j3, the 
reflection coefficient has a minimum in the wave- 
length range between 10 and 30 km. For example, 
when 0=1.0, the optimum wavelength is about 17 
km or approximately 18 kc/s. 

The curves for a 2 in figure 7b also have a particular 
significance. Noting that 

arg R v = — 7r J ra 2 C, 

it is apparent that a 2 C is the phase shift resulting 
from the imperfect reflecting properties of the 
exponential layer. To attach a physical meaning 
to this term it is often desirable to imagine the 
reflection taking place at a height Az x below the 
level z=0. In this case, Az x is chosen so that the 
arg of R v is always — ir. Clearly, 



or 



a 2 = 2kAzi=4TrAz 1 /\ , 



Az l = a 2 \ /(4:T). 



For example, at \ =15 km (i.e., 20 kc/s) and 13=1, 
the effective height of reflection is depressed by 
approximately 3 km. For smaller values of /§, 
it is seen from figure 7b that Az x may be much 
greater. 

6. Final Remarks 

The results, for the coefficients an and a 2 , given 
here may be introduced into the waveguide mode 
theory. In this way, attenuation rates and phase 
velocities of the modes may be obtained. This is a 
valid subject for another paper. 

The influence of nonexponential profiles also 
appears to be a subject worthy of some attention. 
Using the iterative method described in this paper, 
it appears that any reasonably smooth profile 
may be treated in a straightforward manner. Also, 
with some modification, the technique may be used 
to evaluate the effect of discontinuities in gradients 
and local regions of excess ionization. 

General conclusions about the influence of the 
profile of the lower ionosphere on VLF propagation 
must await completion of extensive and systematic 
calculations. It is hoped to report such results 
in the near future. Detailed comparisons of theory 
and experiment are also deferred until these more 
complete computations become available. 
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